REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  Information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington.  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
with  a  collection  of  information  if  it  does  not  display  a  currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. _ _ 

1.  REPORT  DATE  fDD-MM-yYYV9  2.  REPORT  TYPE  3.  DATES  COVERED  fFrom  -  Toj 

07-01-2004  Final  Report  16  August  2002  -  16-Aug-03 


4.  TITLE  AND  SUBTITLE 


Grain-scale  Crystal  Plasticity  Finite  Element  (FE)  Simulations  of  Waspaloy 
During  Hot  Deformation 


6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 

FA8655-02-M4060 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


Professor  Dierk  R  Raabe 


5d.  TASK  NUMBER 


5e.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Max-Planck-lnstitut  fur  Eisenforschung 
Max-Planck-Str.  1 
Duesseldorf  40237 
Germany 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME{S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


EOARD 

PSC  802  BOX  14 
FPO  09499-0014 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


11.  SPONSOR/MONlTOR*S  REPORT  NUMBER(S) 

SPC  02-4060 


14.  ABSTRACT 


This  report  results  from  a  contract  tasking  Max-Planck-lnstitut  fur  Eisenforschung  as  follows:  The  project  aims  to  understand  and  predict  the 
mechanical,  crystallographic,  and  topological  behavior  of  gamma  grains  in  Waspaloy  during  plastic  deformation  at  high  temperature.  The 
study  will  be  conducted  by  use  of  a  crystal  plasticity  finite  element  method  that  updates  the  local  crystallographic  and  hardening  state  of  the 
material  via  Integration  of  the  evolution  equations  for  the  crystal  lattice  orientation  and  the  critical  resolved  shear  stress. 


15.  SUBJECT  TERMS 

EOARD,  Modeling  &  Simulation,  Materials,  Metals  &  alloys,  Nickel,  Superalloy 


20040225  096 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

UNCLAS 


b.  ABSTRACT  c.  THIS  PAGE 

UNCLAS  UNCLAS 


17.  LIMITATION  OF 
ABSTRACT 

UL 


18,  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
OF  PAGES  CHARLES  H.  WARD,  Lt  Col,  USAF 
40 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

+44  (0)20  7514  3154 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39-18 


with  Defense  Federal  Acquisition  Regulation  252  227  701^; 
Declaration  of  Technical  Data  Conformity  (Jan  1997),  252,227-7036, 

“The  Contractor,  Max-Planck-Institut  filr  Eisenforschung,  hereby  delcares  tba*  tn  ti 
belief  the  technical  data  delivered  herewith  under  CoSa  No 
contract  complete,  accurate,  and  complies  with  all  requirements  of  the 


the  requirements  in  Federal  Aquisition  Regulation 
52.227-13,  Patent  Rights — ^Acquisition  by  the  Cf.S.  Government  tJun  19iWi 

(A)  “Disclosures  of  all  subject  inventions  as  defmed  in  FAR  52  227-13  have  been 
reported  in  accordance  with  this  clause  ” 

Or, 

inventions  to  declare  as  defined  in  FAR 
02.227-lo,  during  the  performance  of  this  contract.” 

DATE;  //^  mr 

Name  and  Title  of  Authorized  Oflicial:_  fMf. 


Final  Report  (January,  7*"  2004) 

(Block  17  of  DD  1155,  Accounting  and  Appropriation  Data;  F7LEOA22 120600) 

(AA  5723600  292  47DP  614113  6RER36  592  61 102F  678900) 

CONTRACT  AGREEMENT  NO.  FA8655-02-M-4060 
issued  and  administered  by  EUROPEAN  OFFICE  OF  AEROSPACE,  RESEARCH  AND 
DEVELOPMENT,  223/23 1  OLD  MARYLEBONE  RD,  LONDON  NWl  5TH,  BETH  WANN  0207-514 

X4952,  Beth.Wann@london.af.mil 


Final  Report  on  Grain-Scale 
Crystal  Plasticity  Finite  Element 
Simulations  of  Waspaloy  During 
Hot  Deformation 


Professor  Dierk  Raabe 


Max-Planck-Institut  fur  Eisenforschung 
Max-PlanckStr.  1 
4023 7  Diisseldorf 
Germany 
(raabe@nipie.  de) 


Dierk  Raabe;  Crystal  Plasticity  Simulation  of  Waspaloy  Deformation 


Contents 


1.  Introduction . 3 

2.  Technical  background  and  motivation . 3 

3.  Using  texture  components  in  crystal  plasticity  finite  element  simulations . 4 

3.1  Introduction  to  the  method . 4 

3.2  T  exture  components . 7 

3.3  The  crystal  plasticity  constitutive  model . 8 

3.4  Using  texture  components  in  a  crystal  plasticity  finite  element  simulation . 10 

4.  Description  of  the  model  set-up . 12 

5.  Determination  of  the  constitutive  material  parameters  for  the  compression  simulations . 12 

6.  Results  and  discussion . 14 

6.1  Results  for  perfectly  oriented  cube  crystals  and  different  slip  systems  during  compression 

deformation  (<001  >  axis  parallel  to  compression  axis,  no  initial  scatter) . 14 

6.2  Effect  of  orientational  scatter  prior  to  loading  on  in-grain  texture  evolution  and  grain 

fragmentation  (<001  >  axis  parallel  to  compression  axis) . 16 


6.3  Effect  of  orientational  scatter  prior  to  loading  on  the  lateral  misorientatlon  arrangement  and  on 
the  mechanical  properties  (<001>  axis  parallel  to  compression  axis)  using  {111}<110>  slip  systems...  19 

6.4  Effect  of  orientational  scatter  prior  to  loading  on  the  lateral  misorientatlon  arrangement  and  on 
the  mechanical  properties  (<001  >  axis  parallel  to  compression  axis)  using  {111}[110]  and  {110}[110]- 


type  slip  systems . 23 

6.5  Influence  of  friction  on  in-graIn  texture  evolution  and  grain  fragmentation  for  a  cube-oriented 
grain 

Introduction . 26 

Orientational  stability  of  an  ideal  cube  single  crystal . 27 

Orientational  stability  of  a  cube  single  crystal  with  Gaussian  orientation  scatter . 27 


Orientational  stability  of  a  cube  single  crystal  with  Gaussian  orientation  scatter  under  large-friction  conditions  .  29 
Orientational  stability  of  a  cube  single  crystal  with  Gaussian  orientation  scatter  under  surround  friction  conditions 


. 30 

Influence  of  initial  orientation  scatter  and  of  friction  on  the  re-orientation  behavior  of  cube  crystals . 31 

Conclusions  for  the  friction  study . 32 

6.6  Improved  FE  model  set-up  with  finer  mesh  -  simulation  of  Cube  and  90°  inclined  Cube 
compression . 33 

Mesh  design  and  FE  set-up . 33 

Constitutive  data  fit . 34 

Texture  analysis  after  compression  of  a  Cube  orientation  and  a  Goss  orientation . 35 

7.  Conclusions . 38 


Raabe,  HOARD 


-2- 


MPI  Diisseldorf 


Dierk  Raabe;  Crystal  Plasticity  Simulation  of  Waspaloy  Deformation 


1.  Introduction 

This  report  presents  the  results  obtained  from  crystal  plasticity  FEM  simulations  of  Waspaloy 
compression  deformation  experiments.  The  work  was  conducted  under  the  framework  of  a  mini¬ 
program  on  “Grain-scale  Crystal  Plasticity  Finite  Element  Simulations  of  Waspaloy  During  Hot 
Deformation”.  The  project  was  issued  and  administered  by  the  EUROPEAN  OFFICE  OF 
AEROSPACE,  RESEARCH  AND  DEVELOPMENT  under  the  accounting  and  appropriation 
data  F7LEOA22 120600  (AA  5723600  292  47DP  614113  6RER36  592  61 102F  678900)  and 
contract  agreement  number  FA8655-02-M-4060. 


2.  Technical  background  and  motivation 

The  development  of  polycrystalline,  directionally  solidified,  and  even  single  crystal  superalloys 
and  their  introduction  into  commercial  application  as  turbine  blades  is  one  of  the  major 
technological  achievements  in  materials  science  over  the  past  decades  [1-2].  The  mini-program 
on  “Grain-scale  Crystal  Plasticity  Finite  Element  Simulations  of  Waspaloy  During  Hot 
Deformation”  is  a  numerical  study  which  aims  at  understanding  and  predicting  the 
nonuniformity  of  texture  and  strain  during  hot  compression  of  columnar  grain  structures  having 
different  orientations  relative  to  the  forging  direction.  The  modeling  of  deformation-induced 
plastic  instabilities  is  of  substantial  relevance  for  understanding  the  recrystallization  of 
Waspalloy  during  or  after  hot  compression.  Such  knowledge  is  important  for  obtaining 
homogeneously  recrystallized  material  after  forging  of  cast  blocks  which  have  a  coarse  initial 
grain  structure.  A  refined  homogeneous  microstructure  leads  to  much  better  mechanical 
properties  when  compared  to  coarse  grained  material.  The  current  crystal  plasticity  FEM 
simulation  study  is  in  particular  conducted  for  learning  more  about  the  origin  of  nucleation  sites 
for  dynamic  recrystallization  inside  deformed  grains. 

The  first  part  of  the  study,  which  was  delivered  in  April  2003  was  in  the  form  of  an 
intermediate  report  which  contained  three  main  work  packages.  These  were  first,  the  formulation 
of  a  conceptually  flexible  and  microstmcturally  adequate  crystal-based  FEM  model,  which 
allows  one  to  study  crystallographic  effects  on  the  formation  of  in-grain  orientational  gradients 
in  Waspaloy  during  plastic  compression.  The  model  approach  which  turned  out  to  be  most  suited 
for  this  task  is  the  texture  component  crystal  plasticity  finite  element  method  [12].  Second,  I  had 
to  make  considerations  about  the  potential  slip  systems.  Although  12  slip  systems  of  type 
{lll}<li0>  are  commonly  used  in  crystal  plasticity  simulations  for  FCC  and 
crystallographically  related  materials  recent  single  crystal  experiments  have  shown  that 
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Waspaloy  and  other  FCC  materials  also  tend,  at  elevated  temperatures,  to  deform  additionally  on 
non-octahedral  {110}<110>  slip  systems.  Third,  the  constitutive  behavior  had  to  derived  by 
fitting  the  plastic  behavior  of  Waspaloy  to  mechanical  data  measured  by  Dr.  Semiatin. 

The  consideration  of  non-octahedral  together  with  octahedral  slip  systems  into  the  simulations 
and  fitting  procedures  requires  some  comments  because  one  has  commonly  used  only  12  slip 
systems,  {lll}<liO>,  in  corresponding  crystal  plasticity  simulations  of  the  deformation  of  fee 
single  crystals.  However,  experiments  have  shown  that  single  crystal  Waspaloy  also  tends,  at 
elevated  temperatures,  to  deform  also  on  {110}<li0>  slip  systems.  Therefore,  in  the  present 
crystal  plasticity  simulations  of  high-temperature  deformation  of  cube  oriented  and  related 
coarse  grained  Waspaloy  microstructures  octahedral  and  non-octahedral  slip  systems  are 
considered.  This  approach  of  considering  also  non-octahedral  slip  systems  is  based  on 
experimental  studies  and  numerical  models  of  the  high-temperature  deformation  behavior  of 
single  crystalline  superalloys  [1-4].  For  nickel-base  superalloys  in  general,  and  FCC  single 
crystal  materials  in  particular,  experiments  show  different  operating  slip  systems  such  as 
{111}<110>,  {111}<112>,  and  {110}<110>  [1-2].  Knowles  et  al.  [1]  considered  {111}<112> 
slip  systems  for  a  study  on  creep  of  single  crystal  superalloy  CMSX-4.  Ghosh  et  al.  [2] 
considered  {001  }<1 10>  slip  systems  for  modeling  crystallographic  anisotropy. 


3.  Using  texture  components  in  crystal  plasticity  finite 
element  simulations 

3.1  Introduction  to  the  method 

Metals  mostly  occur  in  polycrystalline  form  where  each  grain  has  a  different  crystallographic 
orientation,  shape,  and  volume  fraction.  The  distribution  of  the  grain  orientations  is  referred  to  as 
crystallographic  texture.  The  discrete  nature  of  crystallographic  slip  along  certain  lattice 
directions  on  preferred  crystallographic  planes  entails  an  anisotropic  plastic  response  of  such 
samples  under  mechanical  loads.  While  the  elastic-plastic  deformation  of  single  crystals  and 
bicrystals  can  nowadays  be  well  predicted,  plasticity  of  polycrystalline  matter  is  less  well 
understood.  This  is  essentially  due  to  the  intricate  interaction  of  the  grains  during  co¬ 
deformation.  This  interaction  leads  to  strong  in-grain  and  grain-to-grain  heterogeneity  in  terms  of 
strain,  stress,  and  crystal  orientation.  One  major  aim  of  polycrystal  research  consequently  lies  in 
identifying  adequate  measures  for  mapping  initial  crystallographic  anisotropy  into  mathematical 
methods  for  predicting  large  strain  plastic  deformation  of  polycrystalline  matter.  The  second 
even  more  demanding  goal  is  the  prediction  of  the  evolution  of  crystalline  anisotropy  during 
plastic  deformation.  This  is  necessary  since  the  crystals  individually  rotate  during  deformation 
owing  to  the  nonsymmetry  of  the  displacement  gradients  created  by  crystal  slip.  The  complex 
microstructural  processes  involved  during  these  reorientation  processes  of  the  individual  grain 
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portions  in  polycrystalline  matter  cannot  be  captured  in  terms  of  simple  empirical  constitutive 
laws  but  require  the  use  of  physically-based  concepts.  In  this  context  the  crystal  plasticity  finite 
element  constitutive  methods  have  increasingly  gained  momentum.  In  these  approaches  one 
typically  assumes  the  stress  response  at  each  macroscopic  continuum  material  point  to  be 
potentially  given  by  one  crystal  or  by  a  volume-averaged  response  of  a  set  of  grains  comprising 
the  respective  material  point.  The  latter  method  naturally  involves  local  homogenization. 
Compared  to  isotropic  J2  approaches  the  crystal  plasticity  finite  element  method  reduces  the 
degrees  of  freedom  for  the  displacement  field  at  each  integration  point  to  the  crystallographic 
slip  dyads  transformed  according  to  the  local  grain  orientation.  Representing  and  updating  the 
crystallographic  orientation  at  each  Gauss  point  renders  crystallographically  discrete  plasticity 
simulations  powerful  tools  for  investigating  anisotropy  and  the  evolution  of  deformation 
textures.  This  chapter  addresses  the  question  how  crystallographic  textures  can  be  merged  with 
crystal  plasticity  finite  element  constitutive  descriptions  in  a  more  rigorous,  scalable,  and 
efficient  way.  A  special  challenge  in  this  context  lies  in  the  reduction  of  redundant  texture 
information  to  a  level  where  sufficient  details  can  be  recovered  without  loosing  physical 
significance. 

This  rather  complex  problem  of  representing  large  microtexture  sets  in  plasticity  simulations  can 
be  split  into  two  separate  tasks.  The  first  one  is  the  formulation  of  a  solution  method  which  uses 
the  crystallographic  orientation  as  a  state  variable.  This  is  typically  achieved  by  formulating  an 
orientation  dependent  constitutive  law  which  maps  the  requested  physical  anisotropy  at  the 
single  crystal  scale  and  by  embedding  this  formulation  into  a  finite  element  code.  The  numerical 
implementation  then  tackles  the  interaction  of  the  differently  oriented  volume  portions  and 
thereby  predicts  the  integral  response  of  the  sample  under  loads.  Any  such  formulation  requires  a 
discrete  representation  of  the  orientation  distribution  function  or  a  portion  of  it  at  each 
integration  point.  Therefore,  the  second  task  consists  in  feeding  one  single  rotation  matrix 
(crystal  orientation)  directly  on  each  Gauss  point  of  the  finite  element  mesh.  This  amounts  to 
mapping  or  respectively  decomposing  orientation  distributions  in  such  a  way  that  they  can  be 
subsequently  mapped  on  a  mesh  in  a  discrete  manner  thereby  matching  the  initial  overall 
distribution. 

One  method  to  map  large  texture  data  sets  onto  a  discrete  mesh  could  be  to  simply  assign  each 
prescribed  or  measured  discrete  orientation  matrix  to  an  individual  grid  point  using.  This  is 
indeed  the  classical  approach  which  is  currently  used  by  many  crystal  plasticity  finite  element 
modelers.  However,  such  an  orientation-by-orientation  (or  grain-by-grain)  approach  is  not 
practicable  when  conducting  forming  simulations  of  larger  parts  typically  containing  more  than 
lO'®  separate  crystals.  The  second  idea,  therefore,  is  to  formulate  a  scalable  method  for  including 
textures  in  crystal  plasticity  finite  element  simulations.  Scalability  means  in  this  context  to  feed 
large  as  well  as  small  texture  data  sets  into  finite  element  models,  but  with  exact  quantification 
of  the  crystallographic  error  arising  from  such  data  reduction. 

Before  discussing  the  advantages  and  disadvantages  of  the  various  functions  available  for  the 
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reproduction  of  textures  it  must  be  noted  that  crystal  plasticity  finite  element  approaches  require 
a  discrete  representation  of  the  orientation  distribution  function  or  a  portion  of  it  at  each 
integration  point.  Mapping  a  discrete  part  of  a  global  texture  to  a  Gauss  point  or  to  a  set  of  Gauss 
points  requires  a  method  which  is  physically  meaningful  (avoiding  over-simplification)  and 
which  at  the  same  time  reduces  the  texture  information  content  to  a  level  which  is  necessary  to 
treat  complex  forming  situations  at  reasonable  computation  costs  (avoiding  over-modeling). 
Such  a  form  is  offered  in  particular  by  the  texture  component  method.  It  approximates  the 
orientation  distribution  function  by  a  superposition  of  sets  of  simple  standard  functions  with 
individual  spherical  coordinates,  orientation  density,  and  scatter  in  orientation  space.  Such  a 
representation  of  a  preferred  orientation  is  referred  to  as  a  texture  component.  In  contrast  to  the 
use  of  global  symmetric  Wigner  functions  for  instance  in  the  Fourier-t^pe  series  expansion 
methods,  the  texture  component  method  is  based  on  using  localized  spherical  normalized 
standard  functions. 

A  new  concept  for  representing  textures  efficiently  in  plasticity  simulations  is  based  on  mapping 
small  sets  of  mathematically  compact  spherical  Gaussian  texture  components  on  the 
integration  points  of  a  crystal  plasticity  finite  element  model.  This  is  done  in  two  steps:  First, 
by  extracting  the  texture  components  and  the  random  background  from  experimental  starting 
textures.  Second,  by  decomposing  these  functions  into  sets  of  discrete  orientations  which 
altogether  exactly  represent  the  original  Gauss  shape  of  each  component  (in  orientation  space)  so 
that  they  can  be  mapped  as  lateral  distributions  onto  a  discrete  lattice.  The  initial  spherical 
distribution  is  thereby  transformed  into  a  spherical  and  lateral  distribution.  During  the 
subsequent  crystal  plasticity  finite  element  simulation  the  texture  component  method  loses  its 
significance  since  each  individual  orientation  originally  pertaining  to  one  of  the  texture 
components  can  undergo  individual  re-orientation  as  in  any  of  the  above  quoted  crystal  plasticity 
methods.  In  order  to  avoid  confusion  one  should,  therefore,  underline  that  the  texture  component 
method  is  used  to  feed  textures  into  finite  element  simulations  on  a  strict  physical  and 
quantitative  basis.  The  components  as  such,  however,  are  not  further  tracked  as  functions  during 
the  simulation.  On  the  other  hand  the  mapped  orientation  points  which  were  extracted  from  the 
components  must  not  be  confused  with  individual  grains,  but  they  are  points  of  a  distribution 
function. 

The  texture  component  approach  ensures  that  the  mapping  of  textures  can  be  linked  to 
materials  concepts  since  texture  components  can  be  directly  interpreted  in  terms  of  characteristic 
metallurgical  mechanisms.  Furthermore  it  makes  sure  that  redundant  texture  information  which 
is  less  relevant,  for  instance  to  predict  plastic  anisotropy,  can  be  reduced  in  the  input  data  set. 
The  following  chapters  explain  the  method  in  more  detail  and  give  examples.  An  overview  of  the 
new  method  is  given  in  [12]. 
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3.2  Texture  components 

The  texture  component  method  offers  a  useful  approach  in  the  context  described  above.  It  is  a 
technique  of  approximating  the  orientation  distribution  function  in  the  form  of  discrete  sets  of 
symmetrical  spherical  model  functions  which  are  defined  in  orientation  space.  The  model 
functions  have  individual  height  and  individual  full  width  at  half  maximum  as  a  measure  for  the 
strength  and  scatter  of  the  crystallographic  texture  component  they  represent.  They  are  typically 
formulated  as  central  functions  which  have  an  isotropic  scatter  in  orientation  space.  The 
mathematical  reproduction  of  the  orientation  distribution  function  by  texture  component 
functions  can  be  expressed  by  the  superposition 

m  =  F  +  'ZPng)  =  j;,rf\g)  where  (1) 

c=l  c=0 

where  g  is  the  orientation,  /(g)  is  the  orientation  distribution  function,  and  F  is  the  volume 
portion  of  all  randomly  oriented  crystals  (random  texture  component).  F  may  be  understood  as 
the  intensity  of  the  only  global  component  used  in  the  approximation,  equivalent  to  /“"(g)  =  1 
for  each  orientation  point  in  Euler  space,  g  e  G .  The  intensity  P  describes  the  volume  fraction 
of  all  crystallites  belonging  to  the  component  c.  The  orientation  density  of  the  component  is 
described  by  a  central  function,  i.e.  its  value  decreases  isotropically  with  increasing  orientation 
distance  =a{g\g)  from  the  maximum.  This  means  that  /'(g)  only  depends  on 
S'  =  ®(g',g)  >  but  it  is  independent  on  the  rotation  axis  n ' . 

The  orientation  distribution  function  is  defined  by 

/ (g)  dg  =  8  which  implies  / (g)  >  0 

where  V  is  the  sample  volume  and  dF^  the  volume  of  all  crystals  with  an  orientation  g  within  the 
orientation  portion  dg=sin((t))  d(j)  d(pi  d(p2.  Normalization  requires 

c 

I /'(g)  dg  =  1  which  implies  ^/'  =  1  (3) 

c=0 

As  a  rule  texture  components  require  positivity,  i.e. 

/'(g)>0  forall  geG  and  /'>0  (4) 

where  G  is  the  orientation  space. 

Spherical  central  functions,  including  corresponding  pole  figures,  can  generally  be  represented  in 
the  form  of  series  expansions  of  x  functions  or  respectively  Legendre  polynomials.  More 
practical  approximations  of  texture  components  have  been  introduced  on  the  basis  of  spherical 
Gauss-  and  Lorentz-fimctions.  The  work  presented  in  this  study  makes  use  of  Gauss-shaped 
model  functions  for  the  decomposition  of  the  orientation  distribution  function  which  are 
described  by 
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P(.S)~  exp(5'  cos  5) 

where 


(5) 


In  2 


and 


(6) 


l-cos(AV2)  /o(5')-/,(S') 

The  according  pole  figure  projections  calculated  in  closed  analytical  form 


^,''(g'>%T)  =  '^'exp(5'sin(uV2))  /o(5'cos(i;V2))  (7) 

where  o'  describes  the  geometry  of  the  component  in  the  respective  pole  figure  projection  and 
I  I  (x)  are  generalized  Bessel  functions.  The  value  is  the  halfwidth  and  can  be  interpreted  as 
the  mean  diameter  of  a  spherical  component  in  orientation  space  and  g°  is  the  center  orientation 
of  the  texture  component. 


The  components  describing can  be  determined  by  the  best  fit  of  the  experimental  pole  figure 
input  data  with  the  recalculated  pole  figures  ^rP^{g\b\y^.  The  index  r 

c 

marks  the  measured  sample  directions  .  The  component  parameters  r ,  g""  and  and  the 
normalization  of  the  pole  figures  are  obtained  by  solving  the  least  squares  problem 


=>  Min. 


(8) 


where  are  weight  factors.  Usually  the  parameters  and  must  be  calculated  by  a  non¬ 
linear  algorithm.  First  estimates  are  required,  which  may  be  obtained  manually  from  the 
graphical  representation  of  the  difference  pole  figures  which  are  calculated  according  to 


A., (>,)=?., («) 

c 

Depending  on  experience  in  interpreting  crystallographic  textures  the  user  can  specify  the 
position,  height,  and  scatter  of  the  texture  components  within  certain  bounds  before  the 
minimization.  This  makes  particularly  good  sense  when  the  number  of  texture  components 
initially  prescribed  to  match  an  experimental  texture  is  small  or  when  a  certain  scatter  width  of 
the  components  should  not  be  exceeded. 


3.3  The  crystal  plasticity  constitutive  model 

In  this  approach  I  use  the  large-strain  constitutive  crystal  plasticity  model  originally  suggested  by 
Kalidindi  [5-8].  In  this  formulation  one  assumes  the  stress  response  at  each  macroscopic 
continuum  material  point  to  be  potentially  given  by  one  crystal  or  by  a  volume-averaged 
response  of  a  set  of  grains  comprising  the  respective  material  point.  The  latter  assumption  can  be 
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referred  to  as  a  local  Taylor-type  or  local  strain-rate  homogenization  assumption.  In  case  of  a 
multi-grain  description  the  volume  averaged  stress  amounts  to 

(T)=.i(w.T.)  00) 

A:=l 

where  N  is  the  total  number  of  individual  orientations  mapped  onto  an  integration  point  using  the 
Taylor  assumption,  the  volume  fraction  of  each  single  orientation  extracted  from  a  texture 
component  as  described  above,  T*  the  Cauchy  stress  produced  by  the  kth  individual  orientation, 
and  (t)  the  volume  average  stress  produced  by  all  orientation  mapped  at  the  integration  point. 
The  constitutive  equation  for  the  stress  in  each  grain  is  then  expressed  in  terms  of 

T*  =  CE*  (*^) 

where  C  is  the  fourth  order  elastic  tensor  and  E*  an  elastic  strain  measure, 

E*  =  -fF*’^  F*-l) 
obtained  by  the  polar  decomposition, 

F  =  F*F'’ 

which  leads  to  a  stress  measure  which  is  the  elastic  work  conjugate  to  the  strain  measure  E*, 

T*  =  F  *■*  (  det(F  *)  T  )  (F 

where  T  is  the  symmetric  Cauchy  stress  tensor  in  the  grain,  and  F*  is  a  local  elastic  deformation 
gradient  defined  in  terms  of  the  local  total  deformation  gradient  F  and  the  local  plastic 
deformation  gradient  F’’.  According  to  equation  (12)  the  relation  between  the  elastic  and  the 
plastic  portion  of  F  amounts  to 

F*  =  F(F'’)r‘  ,  det(F*)>0  ,  det(F'’)=l 

The  plastic  deformation  gradient  is  given  by  the  flow  rule 


with  its  crystalline  portion 


where  m*  are  the  Ath  dyadic  slip  products  of  unit  vectors  b^  in  the  slip  direction  and  normal 
to  the  slip  plane,  and  the  shear  rates  on  these  systems.  The  specific  constitutive  functions  for 
the  plastic  shearing  rates  y^  on  the  slip  systems  are  taken  as 


Yk^Yo 


sgn(rj 


where  is  the  resolved  shear  stress  for  the  slip  system  and  is  the  actual  critical  shear 
stress  on  the  Ath  slip  system.  m  are  material  parameters  representing  shearing  rate  and  the 
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rate  sensitivity  of  slip.  The  calculation  of  has  been  achieved  by  accounting  for  latent 
hardening  through  the  use  of  an  appropriate  hardening  matrix, 

(19) 

i 

where  is  the  rate  of  strain  hardening  on  Ath  slip  system  due  to  a  shearing  on  ith  slip  system, 
is  the  hardening  matrix  describing  the  latent  hardening  behavior  of  a  crystallite,  and  is 
the  hardening  rate  of  the  single  slip  system  i. 


3.4  Using  texture  components  in  a  crystal  plasticity  finite  element 
simulation 

The  main  task  of  the  new  concept  is  to  represent  sets  of  spherical  Gaussian  texture  components 
on  the  integration  points  of  a  finite  element  mesh  for  a  crystal  plasticity  simulation.  This 
procedure  works  in  two  steps:  First,  the  discrete  preferred  orientation  g^  (center  orientation, 
mean  orientation)  is  extracted  from  each  of  the  texture  components  and  assigned  in  terms  of  its 
respective  Euler  triple  ((pi,  (|),  (P2),  i.e.  in  the  form  of  a  single  rotation  matrix,  onto  each 
integration  point  (Fig.  la).  In  the  second  step,  these  discrete  orientations  are  re-oriented  in  such  a 
fashion  that  their  resulting  overall  distribution  reproduces  the  texture  function  which  was 
originally  prescribed  in  the  form  of  a  Gaussian  texture  component.  In  other  words  the  orientation 
scatter  described  initially  by  a  texture  component  function  is  in  the  finite  element  mesh 
represented  by  a  systematically  re-oriented  set  of  orientations,  each  assigned  to  one  integration 
point,  which  reproduces  the  original  spherical  scatter  prescribed  by  that  component.  This  means 
that  the  scatter  which  was  originally  only  given  in  orientation  space  is  now  represented  by  a 
distribution  both,  in  real  space  and  in  orientation  space,  i.e.  the  initial  spherical  distribution  is 
transformed  into  a  spherical  and  lateral  distribution.  The  described  allocation  and  re-orientation 
procedure  is  formulated  as  a  weighted  sampling  Monte  Carlo  integration  scheme  in  Euler  space. 
It  is  important  in  this  context,  that  the  use  of  the  Taylor  assumption  locally  allows  one  to  map 
more  than  one  preferred  crystallographic  orientation  on  each  integration  point  and  to  assign  to 
each  of  them  an  individual  volume  fraction.  This  means  that  the  procedure  of  mapping  and 
rotating  single  orientations  in  accord  with  the  initial  texture  component  scatter  width  is 
individually  conducted  for  all  prescribed  components  as  well  as  for  the  random  background 
extracted  from  initial  experimental  or  theoretical  data. 
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Figure  1  The  decomposition  of  a  cube-oriented  spherical  component  consists  in  extracting  the 
preferred  orientation  from  the  function  (i.e.  (pi=0°,  <|)=0°,  92=0^^  for  cube,  Bunge-Euler  notation) 
and  assigning  it  in  the  form  of  a  single  rotation  matrix  onto  each  integration  point  (left  hand 
side).  In  the  second  step  (right  hand  side)  all  orientation  portions,  defined  by  the  integration 
points,  are  re-oriented  to  reproduce  the  prescribed  orientation  distribution. 


After  decomposing  and  representing  the  initial  texture  components  as  a  lateral  and  spherical 
single  orientation  distribution  in  the  mesh,  the  texture  component  concept  is  no  longer  required 
in  the  further  procedure.  This  is  due  to  the  fact  that  during  the  subsequent  crystal  plasticity  finite 
element  simulation  each  individual  orientation  originally  pertaining  to  one  of  the  texture 
components  can  undergo  an  individual  orientation  change  as  in  the  conventional  crystal 
plasticity  methods.  This  means  that  the  texture  component  method  loses  its  significance  during 
the  simulation.  In  order  to  avoid  confusion  one  should,  therefore,  underline  that  the  texture 
component  method  is  used  to  feed  textures  of  sets  of  sub-grains  or  grains  into  finite  element 
simulations  on  a  strict  physical  and  quantitative  basis.  The  components  as  such,  however,  are  in 
their  original  form  as  compact  functions  not  tracked  during  the  simulation.  It  must  also  be  noted 
that  the  orientation  points  which  were  originally  obtained  fi’om  the  components  do  not  represent 
individual  grains  but  portions  of  an  orientation  distribution  function.  Details  of  this  recently 
introduced  texture  component  crystal  plasticity  finite  element  method  are  given  in  [9-14]. 

The  constitutive  model  described  above  and  the  associated  implicit  time-integration  procedure 
were  implemented  into  the  commercial  finite  element  code  MARC  by  means  of  the  user  defined 
material  subroutine  HYPELA2  [15]. 
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4.  Description  of  the  model  set-up 

The  meshed  sample  was  a  cylinder  with  10  mm  in  diameter  and  15  mm  in  height.  Different 
initial  crystallographic  orientations  parallel  to  the  loading  axis  were  used  for  the  simulations. 
Also,  Compression  was  conducted  at  a  simulated  temperature  of  1450  K  (1177  C),  The  strain 
rate  was  0.005  s\  The  reduction  in  height  amounted  to  50  %  corresponding  to  a  true  strain  of 
about  0.69.  The  friction  in  the  two  contact  planes  were  assumed  at  first  to  be  zero.  A  later  part  of 
the  study  investigates  the  influence  of  changes  in  the  friction  coefficients.  A  set  of  40  eight-node 
3D  brick  elements  was  used  to  discretize  the  mesh  for  the  sample.  The  three-dimensional  finite 
element  mesh  is  shown  in  Fig.  2.  The  compression  tool  was  modeled  as  a  rigid  surface. 


Fig.  2.  The  mesh  used  for  the  finite  element  discretization  of  the  Waspaloy  sample. 
Compression  is  conducted  along  the  z-axis. 


5.  Determination  of  the  constitutive  materiai  parameters 
for  the  compression  simuiations 

The  values  of  the  elastic  parameters  were  taken  to  be  those  of  pure  single-crystalline  nickel  at 
room  temperature,  i.e.  Cii=248.1  GPa,  Ci2=154.9  GPa,  and  C44=124.4  GPa  [16]. 

There  are  two  flow  parameters  in  the  crystal  plasticity  constitutive  equations,  one  is  the  reference 
value  of  slip  rate,  which  is  taken  to  be  0.005  s'',  the  other  is  slip-rate  sensitivity  parameter  m, 
which  is  determined  by  the  macroscopic  strain-rate  jump  experiment  in  simple  compression. 
According  to  the  experimental  data  a  value  of  m  =0.2  was  used  for  the  simulations. 
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The  slip  system  hardening  parameters  were  approximated  by  use  of  an  iterative  fashion  by 
adjusting  the  values  of  these  parameters  in  the  finite  element  simulation  of  simple  compression 
in  such  a  way  that  the  resulting  stress-strain  curve  was  in  good  correspondence  with  the 
experimental  data.  For  12  slip  systems  {lll}<liO>,  the  estimated  initial  value  of  slip  resistance 
So  is  23.8798  MPa  and  the  hardening  parameter  ho  is  0  MPa.  For  18  slip  systems  including  12 
{lll}<liO>  and  6  {110}<li0>  the  estimated  initial  value  of  slip  resistance  so  is  25.6698  MPa 
and  the  hardening  parameter  ho  is  0  MPa.  Fig.  3  shows  the  experimental  data  for  the  stress-strain 
response  and  the  FEM  model  fit  for  that  data. 


Fig.  3.  Experimental  data  for  the  stress-strain  response  and  the  FEM  modePs  fits  for  that  data. 
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6.  Results  and  discussion 


6.1  Results  for  perfectly  oriented  cube  crystals  and  different  slip 
systems  during  compression  deformation  (<001  >  axis  parallel  to 
compression  axis,  no  initial  scatter) 

Figs.  4  and  5  show  the  pole  figures  and  misorientation  contour  plots  for  the  cube  oriented  single¬ 
crystal  Waspaloy  samples  after  deformation.  The  samples  had  no  initial  in-grain  misorientation 
fluctuation,  i.e.  they  were  perfectly  <00 1>  oriented  parallel  to  the  compression  direction  before 
the  compression  deformation.  Fig.  4a  shows  the  pole  figure  results  obtained  for  a  crystal 
plasticity  simulation  with  {111}<110>  slip  systems  (assumption  of  low  temperature 
deformation).  Fig.  4b  shows  the  pole  figures  obtained  for  a  simulation  with  {111}<110>  and 
{110}<110>  slip  systems  (assumption  of  high  temperature  deformation).  The  simulations  shown 
in  Figs.  4  and  5  were  conducted  for  initially  perfect  cube-oriented  single-crystalline  Waspaloy. 
This  means  that  the  texture  component  finite  element  simulations  were  conducted  under  the 
assumption  of  an  orientational  scatter  of  0°  about  the  exact  cube  orientation  prior  to  loading. 
Figs.  4  and  5  show  that  there  is  a  little  lattice  rotation  for  the  {111}<110>  slip  systems  after 
compressing  deformation,  the  maximum  misorientation  is  1.428®C;  but  for  {111}<110>  and 
{ 1 10}<1 10>  slip  systems,  there  is  a  large  lattice  rotation  after  compressing  deformation. 


Single 

orientation 

plots 


001  111  001  111 


iso-intensity 
pole  figure 
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1  1  1 
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RD 


(a)  {1 1 1}<1 10>  slip  systems 


(b)  {111}<110>  and  {110}<110>  slip 
systems 


Fig.  4.  Effect  of  slip  system  selection  on  pole  figures  for  the  cube  orientation  single-crystal 
Waspaloy,  (a)  {1 1 1}<1 10>  slip  systems,  (b)  {111}<110>  and  {110}<110>  slip  systems.  RD 
and  TD  are  given  as  an  arbitrary  in-plane  reference  coordinate  system.  Compression  (50% 
engineering  thickness  reduction)  was  conducted  along  the  normal  direction  ND. 
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(a)  { 1 1 1  }<1 10>  slip  systems  (b)  { 1 1 1  }<1 10>  and  { 1 10}<1 10>  slip  systems 

Fig.  5.  Misorientation  contours  with  (a)  {111}<110>  slip  systems  and  {111}<110>  and 
{1 10}<1 10>  slip  systems,  initial  Cube  orientation. 


Fig.  6  shows  the  contour  plots  of  the  equivalent  total  strains  and  equivalent  Von  Mises  stresses 
for  the  simulations  obtained  for  the  two  different  sets  of  potentially  active  slip  systems.  Fig.  6a 
shows  the  results  obtained  for  {111  }<1 10>  slip  and  Fig.  6b  shows  the  data  for  (111  }<1 10>  and 
{1 10}<1 10>  slip  systems,  respectively.  The  data  show  that  the  distribution  of  the  equivalent  total 
strain  is  similar  for  both  cases.  The  equivalent  Von  Mises  stresses  are  slightly  different  between 
the  two  cases. 


Equivalent  total  strain  Equivalent  total  strain 

(a)  { 1 1 1  }<1 10>  slip  systems  (b)  (111 }<1 10>  and  { 1 10}<1 10>  slip  systems 
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(a)  {111}<110>  slip  systems  (b)  {111}<110>  and  {110}<110>  slip  systems 

Fig.  6.  The  contours  of  equivalent  total  strains  and  equivalent  Von  Mises  stresses  with  (a) 
{1 1 1}<1 10>  slip  systems  and  (b)  {1 1 1}<1 10>  and  {1 10}<1 10>  slip  systems. 


6.2  Effect  of  orientational  scatter  prior  to  loading  on  in-grain  texture 
evolution  and  grain  fragmentation  (<001  >  axis  parallel  to 
compression  axis) 

Figs.  7  and  8  show  the  effects  of  the  initial  orientational  scatter  width  on  the  in-grain  texture 
evolution  for  close-to-cube  oriented  grains  in  terms  of  pole  figures  of  the  deformed  (compressed) 
samples.  The  simulations  were  conducted  for  two  sets  of  potentially  active  slip  systems,  namely, 
for  {111}<110>  slip  systems  (Figs.  7a,  8a)  and  for  {111}<110>  and  {110}<110>  slip  systems 
(Figs.  7b,  8b).  The  pole  figures  show  the  texture  evolution  after  50%  compressing  deformation 
for  the  two  cases. 
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Fig.  7.  Effect  of  the  initial  orientational  scatter  width  on  the  in-grain  texture  evolution 
presented  in  terms  of  pole  figures  after  deformation  with  {111}<110>  slip  systems  (a); 
{111}<110>  and  {110}<110>  slip  systems  (b).  RD  and  TD  are  given  as  an  arbitrary  in-plane 
reference  coordinate  system.  Compression  (50%  engineering  thickness  reduction)  was 
conducted  along  the  normal  direction  ND. 
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Initial 

spherical  Initial  pole  figures  (iso-line  projections) 

scatter  width  of  all  crystals  prior  to  loading. 


Pole  figures  (iso-line  projections) 
after  50%  compressing  deformation 
using  { 11 1  }<1 10>  slip  systems  only. 


Pole  figures  (iso-line  projections)  after 
50%  compressing  deformation  using 
{lll}<110>and  {110}<110>  slip 
systems. 


Fig.  8.  Effect  of  the  initial  orientational  scatter  width  on  the  in-grain  texture  evolution  presented 
in  terms  of  iso-line  pole  figures  after  deformation  with  {111}<110>  slip  systems  (a); 
{111}<110>  and  {110}<110>  slip  systems  (b).  RD  and  TD  are  given  as  an  arbitrary  in-plane 
reference  coordinate  system.  Compression  (50%  engineering  thickness  reduction)  was 
conducted  along  the  normal  direction  ND. 


A  direct  comparison  of  the  pole  figure  data  is  given  in  Fig.  9.  It  documents  the  tendency  of  the 
compressed  Cube  oriented  crystals  that  the  orientational  scatter  evolving  in  them  increases  with 
the  initial  orientation  spread  among  neighboring  grain  portions  and  with  an  increasing  number  of 
glide  systems. 
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Fig.  9.  Direct  comparison  of  data  given  in  Figs.  7  and  8  for  {001 }  pole  figures. 


6.3  Effect  of  orientational  scatter  prior  to  loading  on  the  lateral 
misorientation  arrangement  and  on  the  mechanical  properties 
(<001  >  axis  parallel  to  compression  axis)  using  {111}<110>  slip 
systems 

Figs.  10-12  show  the  effects  of  the  initial  spherical  scatter  width  on  the  lateral  arrangement  of  the 
in-grain  texture  evolution  and  mechanical  properties  for  close-to-cube  oriented  grains.  The 
simulations  were  conducted  for  { 1 1 1  }<1 10>  slip  systems.  The  data  show  the  misorientation  data 
after  50%  compressing  deformation. 
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Fig.  10.  Influence  of  the  initial  orientational  scatter  (5°,  10°, 15°)  on  the  lateral  arrangement  of 
the  in-grain  texture  evolution  for  close-to-cuhe  oriented  grains.  The  simulations  were  conducted 
for  {111}<110>  slip  systems.  The  data  show  the  misorientation  data  after  50%  compressing 
deformation. 
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Lateral  in-grain  distribution  of  von  Mises  stress,  Lateral  in-grain  distribution  of  von  Mises  stress, 
5°  orientation  scatter  prior  to  loading  1 0°  orientation  scatter  prior  to  loading 


Lateral  in-grain  distribution  of  von  Mises  stress, 
15°  orientation  scatter  prior  to  loading 


Fig.  11.  Influence  of  the  initial  orientational  scatter  (5°,  10°,  15°)  on  the  lateral  arrangement  of 
the  in-grain  von  Mises  stress  for  close-to-cube  oriented  grains.  The  simulations  were  conducted 
for  {111}<110>  slip  systems.  The  data  show  the  misorientation  data  after  50%  compressing 
deformation. 
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strain,  10°  orientation  scatter  prior  to  loading 
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Fig.  12.  Influence  of  the  initial  orientational  scatter  (5°,  10°,  15°)  on  the  lateral  arrangement  of 
the  equivalent  in-grain  total  strain  for  close-to-cube  oriented  grains.  The  simulations  were 
conducted  for  {111}<110>  slip  systems.  The  data  show  the  misorientation  data  after  50% 
compressing  deformation. 
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6.4  Effect  of  orientational  scatter  prior  to  ioading  on  the  iaterai 
misorientation  arrangement  and  on  the  mechanical  properties 
(<001  >  axis  parallei  to  compression  axis)  using  {111}[110]  and 
{110}[110]-type  siip  systems 

Figs.  13-15  show  the  effects  of  the  initial  spherical  scatter  width  on  the  lateral  arrangement  of  the 
in-grain  texture  evolution  and  mechanical  properties  for  close-to-cube  oriented  grains.  The 
simulations  were  conducted  for  octahedral  and  non-octahedral  slip  systems.  The  data  show 
the  misorientations  and  mechanical  properties  after  50%  compressing  deformation. 


Lateral  in-grain  misorientation  distribution,  Lateral  in-grain  misorientation  distribution, 
5°  orientation  scatter  prior  to  loading  10°  orientation  scatter  prior  to  loading 


Lateral  in-grain  misorientation  distribution, 
15°  orientation  scatter  prior  to  loading 


Fig.  13.  Influence  of  the  initial  orientational  scatter  (5°,  10°,15°)  on  the  lateral  arrangement  of 
the  in-grain  texture  evolution  for  close-to-cube  oriented  grains.  The  simulations  were  conducted 
for  {111}[110]  and  {110}[110]-type  slip  systems.  The  data  show  the  misorientation  data  after 
50%  compressing  deformation. 
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Fig.  14.  Influence  of  the  initial  orientational  scatter  (5°,  10°, 15°)  on  the  lateral  arrangement  of 
the  in-grain  von  Mises  stress  for  close-to-cube  oriented  grains.  The  simulations  were  conducted 
for  {lll}[liO]  and  {110}[li0]-type  slip  systems.  The  data  show  the  misorientation  data  after 
50%  compressing  deformation. 
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Fig.  15.  Influence  of  the  initial  orientational  scatter  (5°,  10°,  15°)  on  the  lateral  arrangement  of 
the  equivalent  in-grain  total  strain  for  close-to-cube  oriented  grains.  The  simulations  were 
conducted  for  {lll}[li0]  and  {110}[li0]-type  slip  systems.  The  data  show  the  misorientation 
data  after  50%  compressing  deformation. 
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6.5  Influence  of  friction  on  in-grain  texture  evolution  and  grain 
fragmentation  for  a  cube-oriented  grain 

Introduction 

This  chapter  presents  a  study  on  the  sensitvity  of  the  simulation  of  in-grain  orientational 
fragmentation  with  respect  to  friction  conditions.  It  is  particularly  the  Cube  orientation  which  in 
experiments  reveals  substantial  differences  in  its  in-grain  texture  evolution  as  a  function  of  the 
friction  conditions.  In  order  to  scan  various  boundary  conditions  the  simulations  are  conducted 
using  plane  strain  conditions.  The  starting  texture  of  the  Cube  crystals  is  described  in  terms  of  a 
spherical  Gaussian.  Orientational  in-grain  heterogeneity  occurring  during  loading  can,  therefore, 
be  investigated  as  a  function  of  friction  and  orientation  scatter. 

Fig.  16  shows  the  assembly  of  the  finite  element  model.  Four  surfaces  were  used  to  describe  the 
boundary  conditions  for  the  simulated  plane  strain  compression  tests.  This  geometry  and 
individual  surface  treatment  was  chosen  to  mimic  macroscopic  boundary  conditions  typical  of  a 
channel  die  experiment.  The  sample  consisted  of  500  three-dimensional  linear  elements  each 
with  8  integration  points.  The  simulations  used  {111}<110>  slip  systems.  Different  friction 
coefficients  (linear  Coulomb  law)  were  used  to  investigate  the  influence  of  shear.  All  simulations 
were  conducted  to  50%  engineering  thickness  reduction. 


Figure  16  Setup  of  the  finite  element  model  (500  3D  linear  elements).  Four  rigid  surfaces  were 
used  to  simulate  the  boundary  conditions  during  plane  strain  compression,  in  particular  different 
friction  conditions  of  the  different  surfaces. 
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Orientational  stability  of  an  ideal  cube  single  crystal 

Fig.  17  shows  the  simulation  results  obtained  for  the  plane  strain  compression  of  an  initially 
exactly  cube  oriented  and  perfectly  homogeneous  single  crystal,  i.e.  all  elements  had  the  same 
initial  orientation  (pi=0°,  (p2=0°  without  any  initial  orientation  scatter).  The  friction 

coefficients  on  the  top  surface  (contact  surface  with  compression  plane)  and  bottom  surface 
(surfaces  1  and  2  in  Fig.  16)  were  /f=0.1.  The  friction  coefficients  on  the  two  side  surfaces  were 
zero.  Fig.  3  shows  that  the  elements  preserved  a  regular  brick-like  shape  after  50%  deformation 
(engineering  strain  measure).  The  gray  scale  color  code  describes  the  von  Mises  stress.  It  is  very 
homogeneous  throughout  the  specimen.  The  {111}  pole  figure  shows  that  orientational  splitting 
occurred  only  about  ND.  The  orientation  scatter  is  very  symmetric. 
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Figure  17  Simulation  result  for  plane  strain  compression  of  a  cube  single  crystal.  An  exact  cube 
orientation  ((pi=0°,  ((>=0°,  92=0°)  without  any  initial  orientation  scatter  was  assigned  to  all 
elements  prior  to  loading.  The  friction  coefficient  on  surfaces  1  and  2  was  /i=0.1.  The  sample 
was  50%  reduced  in  thickness  (engineering  strain  measure). 


Orientational  stability  of  a  cube  single  crystal  with  Gaussian  orientation 
scatter 

As-crystallized  grains  are  never  perfect  but  contain  lattice  defects  prior  to  loading,  i.e.  the 
orientation  distribution  in  single  grains  is  neither  a  delta  function  in  orientation  space  nor 
laterally  constant.  The  microtexture  in  a  real  single  crystal  can,  therefore,  be  described  in  the 
form  of  an  ideal  (or  exact)  orientation  together  with  an  orientational  scatter  around  it.  In  the 
framework  of  the  texture  component  method  introduced  above  I  describe  the  starting  texture  of 
the  cube  single  crystal  in  the  form  of  a  spherical  Gauss  component  with  a  full  width  at  half 
maximum  of  2.5®.  This  means  that  I  map  a  set  of  single  orientations  onto  the  finite  element  mesh 
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which  satisfies  a  spherical  Gaussian  central  function  about  the  ideal  cube  orientation  (Fig.  18). 
The  500  discrete  orientations  were  assigned  in  random  lateral  order  to  the  elements  of  the  finite 
element  mesh. 


RD 


Figure  18  Gauss-type  spherical  approximation  of  the  initial  microtexture  within  a  cube  single 
crystal  consisting  of  500  single  orientations  in  a  stereographic  {111}  pole  figure  projection.  The 
full  width  at  half  maximum  amounts  to  2.5°. 


Fig.  19  shows  the  deformation  results  for  this  set-up.  The  other  simulation  parameters  were 
identical  to  those  used  in  the  previous  section.  In  clear  difference  to  the  results  for  the  perfect 
initial  cube  orientation  which  shows  that  orientation  splitting  occurs  about  ND  and  RD.  Fig.  22a 
shows  the  corresponding  section  from  the  orientation  distribution  function  (ODF).  It  clearly 
reveals  that  the  orientation  spread  after  deformation  is  larger  about  RD  (-22°  for  f(g)=10)  than 
about  ND  (-14°  for  f(g)=10).  It  must  be  noted  in  this  context  that  the  data  which  are  exactly 
located  on  one  of  the  two  axes  correspond  to  pure  RD  or  respectively  ND  rotations.  The  finite 
element  mesh  and  the  stress  distribution  shows  that  slight  shear  distortion  occurred  in  most 
elements. 


Raabe,  HOARD 


-28- 


MPI  Dusseldorf 


Dierk  Raabe;  Crystal  Plasticity  Simulation  of  Waspaloy  Deformation 


V  Mises  stress 


I+e.5fl0e+02 

+5.91?€+(I^ 
45.£25e+02 
+5.5^9«+(l2 
+5.042e+()2 
44.?50«+»2 
44.45!«+D2 
+4.lS?«+02 
£ — +5.?75«+02 
■ — +$.533e+(l2 
—  -  +3.2S2C+02 

- 43.0fl0e+02 

I— L +2.9374+02 


Figure  19  Simulation  results  for  plane  strain  compression  of  a  single  crystal  with  a  2.5° 
spherical  Gaussian  orientation  scatter  about  the  exact  cube  orientation  prior  to  loading.  Friction 
occurs  only  on  surfaces  1  and  2  (top  and  bottom)  with  a  friction  coefficient  of  /f=0. 1 .  See  also 
the  ODF  section  in  Fig.  22. 


Orientational  stability  of  a  cube  single  crystal  with  Gaussian  orientation 
scatter  under  large-friction  conditions 

Fig.  20  shows  the  simulation  results  for  the  same  case  as  presented  in  Fig.  19,  but  the  friction 
coefficient  amounted  to  /f=0.3.  All  other  parameters  were  the  same  as  in  Fig.  19.  Orientation 
splitting  during  forming  clearly  increases  with  larger  friction  on  the  top  and  bottom  surfaces. 
Although  the  re-orientation  behavior  of  most  elements  is  similar  to  that  shown  in  Fig.  19 
additional  orientation  scatter  appears  about  TD.  The  scatter  about  the  ND  axis  (Fig.  22b)  is 
smaller  than  that  for  the  data  shown  in  Figs.  19  and  22a.  It  is  also  important  to  note  that  the 
overall  orientational  scatter  in  Fig.  20  is  more  smeared  out  than  that  for  the  case  with  smaller 
friction  (Fig.  19).  This  is  also  shown  by  the  ODF  which  reveals  a  much  smaller  maximum  in  the 
orientation  density  at  the  exact  cube  component  (Fig.  22b)  when  compared  to  that  for  the  case 
shown  in  Fig.  22a.  Another  difference  to  Fig.  19  is  that  Fig.  20  shows  more  pronounced  element 
distortion  and  a  more  inhomogeneous  stress  distribution. 


Raabe,  BOARD 


-29- 


MPI  Diisseldorf 


Dierk  Raahe;  Crystal  Plasticity  Simulotion  of  Waspaloy  Deformation 


RD 


V  MIses  stress 

■-+5.91?(+D2 

P<-l-S.042e+«2 
pi-+4.75(l«+d2 
^-■H.45?«+02 
i — +4.1«7e+d2 

- +5.875e+()2 

- +d.5?de+(l2 

-  +d.292€+d2 

- +$.0<IO«+(I2 


Figure  20  Simulation  results  for  plane  strain  compression  of  a  single  crystal  with  a  2.5° 
spherical  Gaussian  orientation  scatter  about  the  exact  cube  orientation  prior  to  loading.  Friction 
only  occurs  on  surfaces  1  and  2  (top  and  bottom)  with  a  friction  coefficient  of  /j=03.  See  also 
the  ODF  section  in  Fig.  22. 


Orientational  stability  of  a  cube  single  crystal  with  Gaussian  orientation 
scatter  under  surround  friction  conditions 

Fig.  21  shows  the  simulation  results  obtained  for  plane  strain  compression  of  a  cube  single 
crystal  with  an  initial  Gaussian  orientation  scatter  of  2.5°  assuming  a  friction  coefficient  of  //=0. 1 
on  all  4  longitudinal  surfaces.  The  other  simulation  parameters  were  the  same  as  in  Figs.  19  and 
20.  The  orientational  scatter  shows  a  strong  component  about  TD  together  with  some  weaker  ND 
and  RD  scatter  (Fig.  22c).  Particularly  the  RD  scatter  is  smaller  than  that  for  the  simulation  with 
non-zero  friction  on  only  two  surfaces  (Figs.  19,  22a). 
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Figure  21  Simulation  results  for  plane  strain  compression  of  a  single  crystal  with  a  2.5"^ 
spherical  Gaussian  orientation  scatter  about  the  exact  cube  orientation  prior  to  loading.  Friction 
occurs  on  4  surfaces  with  a  friction  coefficient  of //=0.1.  See  also  the  ODF  section  in  Fig.  22. 


Influence  of  initial  orientation  scatter  and  of  friction  on  the  re-orientation 
behavior  of  cube  crystals 

Comparison  of  the  results  for  an  ideal  initial  cube  texture  (Fig.  17)  and  for  an  initial  Gaussian 
orientation  scatter  (Figs.  18,  19)  for  identical  friction  conditions  reveals  that  the  simulation 
method  is  sufficiently  sensitive  to  account  for  effects  arising  from  fine  changes  in  starting 
texture. 


Figure  22  Quantitative  presentation  of  the  in-grain  orientation  spread  resulting  after  deformation  in  terms 
of  (p2=0°  sections  from  the  orientation  distribution  function  (ODF)  for  the  three  examples  given  in  Figs.  5- 
7  in  terms  of  pole  figures  (ND:  normal  direction;  RD:  rolling  direction).  The  data  which  are  exactly 
located  on  one  of  the  two  axes  correspond  to  pure  RD  or  respectively  ND  rotations. 

a)  2,5°  orientation  scatter  prior  to  loading,  friction  on  top  and  bottom,  //=0. 1 

b)  2.5°  orientation  scatter  prior  to  loading,  friction  on  top  and  bottom,  //=0.3 

c)  2.5°  orientation  scatter  prior  to  loading,  friction  on  4  surfaces,  //=0.1 
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These  results  confirm  many  experiments,  namely,  that  the  cube  orientation  is  metastable  under 
plane  strain  conditions.  This  was  also  found  in  recent  theoretical  studies  on  the  re-orientation 
behavior  of  different  orientations  by  use  of  homogenization  and  crystal  plasticity  finite  element 
analysis  [17].  This  investigation  also  confirmed  that  orientation  gradients  in  initially  uniform 
cube  crystals  can  occur  under  gradient-free  external  loads.  The  intrinsic  origin  of  this  effect  was 
quantified  in  terms  of  the  change  in  crystal  re-orientation  upon  small  changes  in  initial 
orientation.  Such  starting  conditions  are  in  the  present  study  represented  by  the  Gaussian 
orientation  scatter  prior  to  deformation.  The  theory  was  formulated  as  a  divergence  operator 
applied  to  re-orientation  rate  vector  fields.  The  obtained  scalar  divergence  function  gave  an 
excellent  quantification  of  the  stability  of  grains  under  homogeneous  boundary  conditions  as  a 
function  of  their  orientation.  Positive  divergence  which  is  a  source  in  the  re-orientation  rate 
vector  field  characterize  orientations  with  diverging  non-zero  reorientation  rates  which  are 
kinematically  instable  and  prone  to  build  up  orientation  gradients.  Besides  such  basic 
kinematical  interpretation  I  attribute  the  deviation  between  Figs.  17  and  19  also  in  terms  of  the 
loss  in  local  symmetry  (Fig.  19).  The  slight  element  distortion  shows  that  due  to  the  laterally 
random  distribution  of  the  imposed  initial  2.5°  orientation  scatter  the  plane  strain  state  is  locally 
violated  entailing  strong  RD  rotations  (Fig.  22).  These  observations  correspond  to  recent  results 
which  showed  that  weakly  positive  divergent  orientations,  such  as  cube  in  FCC,  reveal  a  larger 
dependence  of  their  re-orientation  behavior  on  the  local  mechanical  boundary  conditions  (local 
deviation  from  a  plane  strain  state)  and,  therefore,  a  larger  variety  of  the  resulting  orientation 
spreads,  than  less  divergent  crystals.  Such  effects  arising  from  the  loss  in  local  symmetry  are 
surely  assisted  by  friction.  However,  the  influence  of  friction  is  likely  to  promote  TD  rotations 
rather  than  RD  rotations  because  it  induces  longitudinal  shear  strains.  With  increasing  friction 
coefficient  and  increasing  number  of  surfaces  imposing  non-zero  friction,  the  orientation  and 
stress  distribution  becomes  increasingly  heterogeneous. 

Conclusions  for  the  friction  study 

This  section  discussed  the  orientational  stability  of  Cube  FCC  crystals  under  plane  strain  (50% 
engineering  thickness  reduction)  by  use  of  a  texture  component  crystal  plasticity  finite  element 
method  (12  slip  systems).  The  main  results  are: 

•  If  the  friction  is  small  (//=0.1)  and  the  initial  orientation  of  the  sample  is  exactly  cube 

everywhere  in  the  mesh  (all  elements  have  the  same  initial  orientation  (pi=0°,  (t)=0°, 
(p2=0°  without  any  initial  orientation  scatter  prior  to  loading)  the  orientation  spread  after 
deformation  can  be  described  in  terms  of  a  pure  ND  rotation. 

•  If  the  friction  is  small  (//=0.1)  and  the  initial  orientation  of  the  sample  is  a  cube  orientation 

with  2.5°  spherical  Gaussian  scatter  the  orientation  spread  after  deformation  can  be 
described  in  terms  of  ND  and  RD  rotations. 

•  If  the  friction  is  large  (/^0.3)  and  the  initial  orientation  of  the  sample  is  a  cube  orientation 

with  2.5°  spherical  Gaussian  scatter  the  orientation  spread  after  deformation  can  be 
described  in  terms  of  ND,  RD,  and  TD  rotations.  The  latter  rotation  mode  is  due  to  the 
induced  forward  shear. 
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6.6  Improved  FE  model  set-up  with  finer  mesh  -  simuiation  of  Cube 
and  90''  inclined  Cube  compression 

Mesh  design  and  FE  set-up 

The  second  approach  for  an  improved  meshing  used  again  a  cylinder  with  10  mm  in  diameter 
and  15  mm  in  height.  Different  initial  crystallographic  orientations  parallel  to  the  loading  axis 
were  used  for  the  simulations.  Compression  was  conducted  at  a  simulated  temperature  of  1450  K 
(1177  C),  The  strain  rate  was  varied  between  0.005  s’^  and  0.1  s'^  The  reduction  in  height 
amounted  to  50  %  corresponding  to  a  true  strain  of  about  0.69.  The  friction  in  he  two  contact 
planes  were  assumed  to  be  zero. 


Fig.  23.  The  improved  mesh  used  for  the  finite  element  discretization  of  the  Waspaloy  sample. 
Compression  is  conducted  along  the  z-axis. 
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Fig.  24.  Constitutive  fit  for  the  Cube  orientation.  Stress-strain  response  and  FEM  model  fits. 


1 


Raabe,  HOARD 


-34- 


MPI  Dusseldorf 


Dierk  Raabe;  Crystal  Plasticity  Simulation  of  Waspaloy  Deformation 


Fig.  25.  Constitutive  fit  for  the  Goss  orientation  (90°  TD  rotated  Cube  orientation).  Stress-strain 
response  and  FEM  model  fits. 


Texture  analysis  after  compression  of  a  Cube  orientation  and  a  Goss 
orientation 

Figs.  26-28  show  the  pole  figures  for  the  Cube  and  the  Goss  oriented  single-crystal  Waspaloy 
samples  after  compression  deformation.  Fig.  4a  shows  the  pole  figure  results  obtained  for  a 
crystal  plasticity  simulation  with  {111}<110>  slip  systems  (assumption  of  low  temperature 
deformation).  Fig.  26  shows  the  pole  figures  after  compression  for  the  Cube  orientation.  The 
center  of  the  pole  figures  corresponds  to  the  compression  direction.  Fig.  27  reveals  the  pole 
figures  after  compression  for  the  Goss  orientation.  The  simulations  were  conducted  using  12  slip 
systems  {111}<110>.  Fig.  28  shows  the  corresponding  data  for  the  Goss  orientation  using  18 
slip  systems  {111}<110>  and  {100}>110>.  The  corresponding  microstructures  are  presented  in 
Figs.  29-31. 
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Fig.  26.  Pole  figures  after  compression.  Center  of  pole  figure  corresponds  to  the  compression 
direction.  Cube  orientation,  same  result  for  12  slip  systems  {111}<110>  as  well  as  for  18  slip 
systems  {111}<110>  and  {100}>110>. 
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Fig.  27.  Pole  figures  after  compression.  Center  of  pole  figure  corresponds  to  the  compression 
direction.  Goss  orientation,  12  slip  systems  {1 1  !}<!  10>. 
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Fig.  28.  Pole  figures  after  compression.  Center  of  pole  figure  corresponds  to  the  compression 
direction.  Goss  orientation,  18  slip  systems  {1 1 1}<1 10>  and  {100}>1 10>. 
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Fig.  29,  Lateral  in-grain  misorientation  distribution, 

Cube  orientation,  same  result  for  12  slip  systems  {111}<110>  as  well  as  for  18  slip  systems 
{111}<110>  and{100}>110>. 
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Fig.  30.  Lateral  in-grain  misorientation  distribution, 
Goss  orientation,  12  slip  systems  { 1 1 1  }<1 10>. 
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Fig.  31.  Lateral  in-grain  misorientation  distribution, 

Goss  orientation,  18  slip  systems  {1 1 1}<1 10>  and  {100}>110>. 


7.  Conclusions 

The  current  study  investigated  the  orientation  dependence  of  the  in-grain  misorientation 
development  of  Waspaloy  as  a  function  of  the  starting  orientation.  It  was  found  that  orientations 
which  are  inclined  with  respect  to  the  Cube  orientation  reveal  a  much  stronger  tendency  to 
undergo  orientational  grain  fragmentation  than  Cube  oriented  grains.  The  Goss  orientation 
reveals  a  strong  tendency  towards  orientational  grain  fragmentation  when  compared  to  the  Cube. 
Additionally  the  influence  of  friction  on  such  phenomena  was  systematically  studied.  It  was 
observed  that  larger  friction  coefficients  lead  to  a  stronger  in-grain  misorientation  development. 
This  aspect  is  of  relevance  particularly  in  cases  where  Cube  oriented  as-cast  grains  cannot  be 
rotated  before  the  deformation. 
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